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^ ■ Abstract 

We study a number of models, based on a non-Abelian discrete group, that 
successfully reproduce the simple and predictive Yukawa textures usually as- 
sociated with U(2) theories of flavor. These models allow for solutions to the 
solar and atmospheric neutrino problems that do not require altering suc- 
cessful predictions for the charged fermions or introducing sterile neutrinos. 
Although Yukawa matrices are hierarchical in the models we consider, the 
| mixing between second- and third-generation neutrinos is naturally large. We 

first present a quantitative analysis of a minimal model proposed in earlier 
work, consisting of a global fit to fermion masses and mixing angles, including 
^ . the most important renormalization group effects. We then propose two new 

variant models: The first reproduces all important features of the SU(5) xU(2) 
unified theory with neither SU(5) nor U(2). The second demonstrates that 
discrete subgroups of SU(2) can be used in constructing viable supersymmet- 
ric theories of flavor without scalar universality even though SU(2) by itself 
cannot. 
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I. INTRODUCTION 



It is possible that the observed hierarchy of fermion masses and mixing angles originates 
from the spontaneous breakdown of a new symmetry Gf that acts horizontally across the 
three standard model generations. Ideally, all Yukawa couplings except that of the top quark 
are forbidden by Gf invariance at high energies; the remaining ones are generated when a set 
of fields that transform nontrivially under Gf develop vacuum expectation values (vevs). 
A hierarchy in couplings is obtained if Gf is broken sequentially at energy scales /ij through 
a series of nested subgroups , such that 

Gf J^ Hl J^H 2 J^... for /i x > // 2 > // 3 • • • • (1.1) 

At each stage of the symmetry breaking there is an associated small dimensionless parameter 
((f)i)/Mf, where 0j is a 'flavon' field whose vev is responsible for the breaking Hi_i — > Hi, 
and where Mf is the ultraviolet cutoff of the G/-invariant effective theory. The ratios <f>i/Mf 
appear in higher- dimension operators that contribute to Yukawa couplings in the low-energy 
theory. For example, the superpotential term 

-LQ 3 H D cf> b D 3 (1.2) 

leads to a bottom quark Yukawa coupling of order (</>&)/M/. The most general set of opera- 
tors involving the fields of the minimal supersymmetric standard model (MSSM) and the <fi 
fields must provide for Yukawa textures that are phenomenologically viable. If flavor univer- 
sality of scalar superpartner masses is not simply a consequence of the mechanism by which 
supersymmetry breaking is mediated fl|-§J, then a successful model must also explain why 
these scalars do not contribute to flavor- changing neutral current processes at unacceptable 
levels. 

Models with horizontal symmetries have been proposed with Gf either gauged or 
global, continuous or discrete, Abelian or non-Abelian, or some appropriate combination 
thereof HH- Abelian flavor symmetries have been used successfully to explain the absence 
of supersymmetric flavor-changing processes by aligning the fermion and sfermion mass ma- 
trices ||. However, the freedom to choose a number of new U(l) charges for each MSSM 
matter field represents so much freedom that these models seem ad hoc, at least from a low- 
energy point of view. Non-Abelian symmetries are more restrictive, as the Yukawa matrices 
generally decompose into a smaller number of irreducible Gf representations. Thus, it is not 
unreasonable to expect that minimal models exist that are both successful and aesthetically 
compelling. This is the primary motivation for the current work. 

In non-Abelian flavor models, the existence of three generations of matter fields, the 
heaviness of the top quark, and the absence of supersymmetric flavor- changing processes 
together suggest a 2©1 representation structure for the MSSM matter fields. With this 
choice it is not only possible to distinguish the third generation, but also to achieve an exact 
degeneracy between superparticles of the first two generations when Gf is unbroken. In the 
low-energy theory, this degeneracy is lifted by the same small symmetry-breaking parameters 
that determine the light fermion Yukawa couplings, so that FCNC effects remain adequately 
suppressed, even with superparticle masses less than a TeV. 
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A particularly elegant model of this type considered in the literature assumes the con- 
tinuous, global symmetry Gf =U(2) |8HlO|. Quarks and leptons are assigned to 2©1 repre- 
sentations, so that in tensor notation one may represent the three generations of any matter 
field by F a + F 3 , where a is a U(2) index, and F is Q, U, D, L, or E. A set of flavons is 
introduced consisting of 4> a , S a b, and A a b, where is a U(2) doublet, and S(A) is a symmetric 
(antisymmetric) U(2) triplet (singlet). The doublet and triplet flavons acquire the vevs 

and #=(°°), (1.3) 
M, \ M, \0 ej' ^ 

the most general set of nonvanishing entries consistent with an unbroken U(l) symmetry 
that rotates all first generation-fields by a phase. This residual U(l) symmetry is broken at 
a somewhat lower scale by the flavon A 



(A) 




:i-4) 



where e' < e. Thus, the sequential breaking 

U(2) -U U(l) nothing, (1.5) 
yields a Yukawa texture for the down quarks, for example, of the form 

d x e' 

Y D « I -d.e' d 2 e d 3 e I £ , (1-6) 



<i 4 e d 
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where di, . . . , d 5 are 0(1) coefficients. With the choice e ~ 0.02 and e' ~ 0.004, this texture 
achieves the correct hierarchy in down quark mass eigenvalues and gives contributions of the 
appropriate size to entries of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The 0(1) 
coefficients may be determined from a global fit, as in Ref. fLOf . The ratio mf,/m t is assumed 



to be unrelated to U(2) symmetry breaking, and is simply put into the low-energy theory 
by hand. This is accomplished by choosing the free parameter £ in Eq. ( |1.6| ). 

While the form of Yd is viable, U(2) symmetry by itself cannot explain the differences 
between the hierarchies within Yd and Yy. Quark mass ratios renormalized at the grand 
unified scale are given approximately by 

rrid m s :: rn& = A 4 :: A 2 :: 1, (1.7) 

while 

m u :: m c :: m t = A 8 :: A 4 :: 1, (1.8) 

where A ~ 0.22 is the Cabibbo angle. Clearly, an additional suppression factor p is required 
in Y[/ for those elements that contribute most significantly to the up and charm quark mass 
eigenvalues, 

/ uxe'p 
Yjj w -me'p u 2 tp u 3 e \ . (1.9) 
\ M 4 e M5 
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where U\ . . . U5 are 0(1) coefficients. By embedding the U(2) model in a grand unified theory- 
it is possible to obtain p « e naturally; the model can then accommodate all the desired 



fermion mass hierarchies for choices of the coefficients Ui and dj that are all of order one flit) 
For example, in an SU(5) GUT, Yy is associated with the coupling 10-10-5, where the 10's 
represent matter fields, and the 5 is the Higgs field H . However, 

10® 10 = 5,©45 a ©50 s , (1.10) 

where the subscripts indicate symmetry or antisymmetry under interchange of the two 10's. 
If we assume that the antisymmetric flavon A is an SU(5) singlet, the product AH is a 
5 tt , and does not contribute to Yy. Similarly, if the flavon S is a 75 with a vev in the 
hypercharge direction in SU(5) space, then the part of SH that contains the Higgs doublet 
field transforms as a 45 s , which again does not contribute to Yy. To obtain nonvanishing 
couplings of the right size in the upper two-by-two block of Yy one introduces a singlet 
flavon E that transforms as an SU(5) adjoint. The vev of S implies that the breakings of 
both U(2) to U(l) and SU(5) to the standard model gauge group are associated with vevs 
of order e. Thus, it is natural to assume (E) ~ e, which provides exactly the desired value of 



p in Eq. (1.9). Moreover, the SU(5) assignments for A and S provide for a Georgi-Jarlskog 



mechanism [|ITJ, so that unified U(2) models successfully account for the charged lepton 
mass spectrum as well. 

While the textures that follow from the simple two-step breaking of a U(2) flavor sym- 
metry are indeed minimal, the original symmetry group is not. It is natural to ask whether 
there are small discrete groups that work equally well as horizontal symmetries. It was 
shown in Ref. fll2[ that the charged fermion Yukawa textures usually associated with U(2) 
models may be reproduced assuming the symmetry Gf — T' x Z 3 , and the breaking pattern 

T' <g> Z 3 — ^ Z° nothing. (1.11) 

Here, T' is the double tetrahedral group, a discrete subgroup of SU(2) corresponding to the 
symmetry of a regular tetrahedron. The factor is the diagonal subgroup of a Z3 C T' 
and the additional Z 3 factor (see Section [TV]). Since U(2) is isomorphic to SU(2)xU(l), it 
is not surprising that our discrete symmetry is a product of a discrete subgroup of SU(2) 
and a discrete subgroup of U(l). Moreover, it was argued in Ref. [jnj that this symmetry is 
minimal in the sense that 

• T' is the smallest discrete subgroup of SU(2) (and in fact the smallest group of 
any kind) with 1-, 2- and 3-dimensional representations and the multiplication rule 
2®2=3©1. These two ingredients are necessary to reproduce the successful U(2) 
textures. 

• Z 3 is the smallest discrete subgroup of U(l) that allows Gf to contain a subgroup 
forbidding all order O(e') entries in the Yukawa textures. 

The latter statement applies to models in which T' is a discrete gauge symmetry (see Sec- 
tion U); models with a global T" symmetry do not require any additional Abelian factors, as 
we demonstrate in Section |V11|. The use of discrete gauge rather than global symmetries is 



motivated by various arguments that the latter are violated at order one by quantum grav- 



itational effects fll^l - In two of the models we present, T' is an anomaly-free discrete gauge 
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symmetry, while the additional Z n factor is not. As in many of the Abelian models described 
in the literature ||, we simply assume that the Z n factor may be embedded in a U(l) gauge 
symmetry whose anomalies are cancelled by the Green-Schwarz mechanism WM. Thus, our 



models may be viewed as consistent low-energy effective theories for flavor symmetries that 
are local in a complete, high-energy theory. 

On a more practical level, the different representation structure of T' allows for elegant 
solutions to the solar and atmospheric neutrino problems that do not alter the predictive 
quark and charged lepton Yukawa textures, nor require the introduction of sterile neutrinos. 
While similar results can be obtained in some SO(10)xU(2) models [|15| , we obtain our 
successful solutions using a much smaller symmetry structure.[] One goal of this paper is to 
study these solutions at a level of quantitative detail not presented in our earlier work. 

In addition, we propose two new models involving T' symmetry. The first model, based 
on the discrete gauge symmetry T'xZq, reproduces all important features of the SU(5) x U(2) 
model without requiring a field-theoretic grand unified theory. In other words, the suppression 
of m u and m c in the SU(5)xU(2) theory described earlier is achieved in T' x Z 6 without 
SU(5). In addition, the ratio mb/m t , which is not explained in SU(5)xU(2), is predicted 
in our model to be of 0(e) w 0.02 for tan/3 ~ 0(1), where tan/3 is the ratio of Higgs 
field vevs (Hjj)/(H d ). In a second model, we consider the implications of T' as a purely 
global flavor symmetry. Although in this case the symmetry may not be fundamental, it 
could still arise as an accidental symmetry at low energies. We show that it is possible to 
construct a viable model based on T' alone, with no additional Abelian factors. While it 
is well known that supersymmetric models with a continuous SU(2) flavor symmetry and 
a 2 © 1 representation structure do not have viable Yukawa textures, our global T' model 
demonstrates that discrete subgroups of SU(2) remain viable alternatives. 

Our paper is organized as follows. In the next section, we discuss the meaning of discrete 
gauge symmetries and the relevant anomaly-cancellation constraints in the low-energy effec- 
tive theory. In Sections [TTT| and |TV|, we review the group theory of T' and present the minimal 
model described in Ref . |12| . In Section [Vj we fit predictions of the model to charged fermion 
and neutrino masses and mixing angles, including the most significant renormalization group 
effects. In Section [VI], we present the T' x Z§ model that reproduces the important features 
of the SU(5) xU(2) model with neither SU(5) nor U(2). In Section |VII| , we show how to con- 



struct a viable global T' model with no Abelian factors. In Section |V111] we comment on one 



scenario involving sterile neutrinos, and in the final section we summarize our conclusions. 



II. WHAT IS A DISCRETE GAUGE SYMMETRY? 

Let us define a discrete gauge symmetry provisionally as a discrete remnant of a sponta- 
neously broken continuous gauge symmetry. Below the breaking scale A of the continuous 
symmetry, the low-energy effective Lagrangian has interactions that are invariant under the 
unbroken discrete group, no massless gauge fields, and derivatives that transform trivially. 
It would seem then that this effective theory is identical to one with a purely global discrete 



1 For a similar approach, see Ref. [16|. 
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symmetry. In this section, we review the arguments suggesting that this is not the case. We 
first illustrate how gauge invariance of a theory spontaneously broken to a discrete subgroup 
dictates the form of all terms in the low-energy effective theory, and thus renders its dis- 
crete invariance immune to wormhole dynamics. We then show that a theory with a discrete 
gauge symmetry predicts topological defects not present in a theory with a global symmetry, 
and that these play an important role in demonstrating that discrete gauge charges leave 
quantum-mechanical hair on black holes. Both observations suggest that discrete gauge 
symmetries are viable as candidates for fundamental symmetries of nature. After reviewing 
these arguments we summarize the anomaly-cancellation constraints relevant to low-energy 
theories with discrete gauge symmetries. We use these constraints in constructing models 
throughout this paper. 

Following a discussion by Banks , let us consider the low-energy effective theory that 
results from spontaneously breaking a U(l) gauge symmetry to a discrete subgroup. The 
full theory consists of two scalar fields x an d 4> with U(l) charges q and 1, respectively, 
where q is an integer. The Lagrangian is the usual one for an Abelian Higgs model: 

C = -^F, U F^ + \d, X ~ nKx\ 2 + \d,4> - iA^\ 2 + V{ X ] X)- (2.1) 

The potential V is such that the x field acquires a vacuum expectation value A. Let us 
rewrite the Lagrangian using the nonlinear field redefinition x = (A + cr)e ,e / \/2. This yields 

C = —F^F^ + + i(A + <t) 2 (3„0 - qA,) 2 + \d^ - iA^\ 2 + V(a) , (2.2) 

where a is the Higgs field and 9 is the would-be longitudinal component of the U(l) gauge 
boson in unitary gauge. We choose to construct a low-energy effective theory in which the 
a field, which has a mass of order A, is integrated out. However, we retain the gauge field 
A M as well as the unphysical scalar field 9. Although the gauge symmetry is spontaneously 
broken, the Lagrangian of the theory remains invariant under the local U(l) transformation: 

0_> e <«(*)^ A» ^> AO + dpOtix), 9^9 + qa(x). (2.3) 

The low-energy effective Lagrangian then consists of the kinetic terms 

C = + - iA ^? + - qA,) 2 , (2.4) 

as well as the most general set of gauge-invariant operators involving the fields 0, e te , and 
covariant derivatives, with powers of A included to obtain the correct mass dimensions. We 
can classify the interactions in the effective Lagrangian that involve into two types: terms 
that are invariant under global U(l) transformations on <fr alone (with the other fields held 
fixed), and those that are not. A typical term of the first type is (jy<j)\ terms of the second 
type necessarily involve the U(l) gauge-invariant product 

e~V , (2.5) 

or similar products with derivatives. Such terms are invariant under a Z q phase rotation 
of the field alone. Thus, gauge invariance of the low-energy theory implies it must have 
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an unbroken Z q symmetry. Since this is a consequence of a local symmetry, it cannot be 
violated by wormhole dynamics. 

We now show that information on discrete gauge charges is not lost when a charged 
particle falls into a black hole. To do so, first note that the Abelian Higgs model has stable 
cosmic string solutions. In the case where = 0, the kinetic energy terms in Eq. ( |2.4| ) are 
minimized when 

A, = -d,9 . (2.6) 
q 

For nonsingular gauge field configurations, this is related to A^ = by a gauge transforma- 
tion. However, singular solutions also exist; a cosmic string along the x 3 axis corresponds 
to 

1 x j 

Ai = -eij—^ k , i, j = 1,2 , 9 = arctanfe/xi) . (2.7) 

q xj + x\ 

If one couples the gauge field to a classical current j^, then the change in the action by 
adding one such cosmic string is 

5S = (1/q) J <W , (2.8) 

which follows from Eq. ( |2.6|) . Taking j M to be the current of a particle with unit U(l) charge 
(and hence nontrivial Z q charge) that circles the string, one finds that 

5S = — . (2.9) 
Q 

This implies an observable Aharanov-Bohm effect in the scattering of particles with discrete 



gauge charge off cosmic strings. Krauss and Wilczek [18] use this observation to argue that 
the scattering of a cosmic string off a particle with discrete gauge charge that is falling into 
a black hole is insensitive to the point at which the particle crosses the event horizon. Thus, 
the discrete charge of the particle is not lost, and the black hole grows quantum- mechanical 
hair. 

It is interesting to note that the discussion above may be rephrased in unitary gauge by 
making the initial replacements 

B li = A li -(l/q)d li 9, and <£> = e~ ie ^cj> , (2.10) 

in Eq. (|2.2j ), which then becomes 

47 F ^ + ^ + ? 



C = -— F ^ + nfaVv + o ( A + °?q 2 B»B, + \d^ - iB^\ 2 + V{a) . (2.11) 



Unlike the previous approach, all the fields above are gauge-invariant; one may integrate 
out B^ and a, and obtain all the possible Z g -invariant interactions involving the light field 
$. This formulation of the low-energy theory is peculiar in that the periodicity of 9 implies 
that 



g2n-7ri 



J /9$ = $ j f or a u integers n. (2.12) 
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Thus, the field manifold of is not the complex plane C, but rather the orbifold C/Z q : Field 
configurations connected by Z q transformations are identified, and hence are physically 
redundant, the hallmark of a gauge symmetry. Given this manifold, the field $ has a 
conical singularity at the origin in field space; strings in unitary gauge correspond to $ field 
configurations that wrap around this singularity as the azimuthal angle varies from to 2n. 

As the previous U(l) — > Z q example demonstrates, a discrete gauge symmetry can arise 
in a renormalizable field theory when a continuous gauge symmetry is spontaneously broken 
by a Higgs field vev that leaves a discrete symmetry unbroken. The same can occur for 
non-Abelian symmetries as well. For example, one may break a gauged SU(2) symmetry 
with a Higgs field transforming as a 7 (which contains a T" singlet), leaving the theory 
invariant under X". On the other hand, the U(l) — > Z q example suggests how a discrete 
symmetry may be defined without an explicit embedding in a continuous group. In string 
theory, the discrete symmetry may be a remnant of general coordinate invariance, ordinary 
gauge invariance, or the larger gauge symmetry of string theory [JT7| . For our purposes, 
however, the nature of the high energy theory is irrelevant. 

It is worth mentioning in passing that spontaneously-broken discrete gauge symmetries 
have domain walls that are not topologically stable. Holes bounded by strings may spon- 
taneously nucleate, allowing the walls to tear themselves to pieces while dissipating energy 
through gravitational radiation [IS]. The effectiveness of this mechanism at avoiding cosmo- 
logical problems is not relevant to our discussion since the flavor-symmetry-breaking scale 
in our models is high enough (of order the unification scale) that all topological defects are 
eliminated by inflation. 

Finally, it is relevant to consider whether there are any constraints on the low-energy 
particle content of theories with discrete gauge symmetries. Since continuous gauge sym- 
metries must satisfy anomaly-cancellation conditions, the particle content of low-energy 
theories with discrete gauge symmetries is restricted. Ibanez and Ross [20| were the first to 
consider the constraints on a discrete gauged Z q symmetry, and their results were refined 



by Banks and Dine [21 1: Let Go be a simple factor of the continuous group in which a dis- 



crete gauge symmetry is embedded, and let Ga and G^ represent the unbroken Abelian and 
non-Abelian gauge symmetries of the low-energy effective theory. Cancellation of the GqG 2 n 
anomaly is the only new requirement for consistency of the low-energy theory; all other 
anomaly-cancellation constraints involving Go can be satisfied by the introduction of heavy 
states. Banks and Dine point out that this requirement, termed the linear Ibanez- Ross con- 
dition, is equivalent to demanding discrete gauge invariance of nonperturbative interactions 
generated by instantons of the unbroken continuous gauge groups. This observation demon- 
strates that consistency of a discrete gauge symmetry at low energies can be established 
without reference to any particular embedding. 
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III. THE GROUP T' 



All of the symmetries described in this paper contain T", the double tetrahedral group. [] 
Geometrically, T' is defined as the group of all 24 proper rotations in three dimensions leaving 
a regular tetrahedron invariant in the SU(2) double covering of SO(3). This perhaps opaque 
definition may be understood in the following way. There exists a group of 12 elements 
called the tetrahedral group T, consisting of all proper rotations in three dimensions leaving 
a regular tetrahedron invariant (Fig. 1). It is constructed by parameterizing the group 
SO (3) of all proper rotations in three dimensions in terms of familiar Euler angles, and then 
restricting to those discrete values of angles describing rotations taking a regular tetrahedron 
into coincidence with itself. The same Euler angles describe rotations in SU(2) space, since 
SU(2) and SO(3) are locally isomorphic, so that T' is the subgroup of SU(2) corresponding 
to the same Euler angles as T C SO (3). One therefore expects that even- dimensional 
representations of T' are spinorial, i.e., are multiplied by —1 under a 2n rotation (called R 
in the literature), while odd- dimensional representations of T' coincide with those of T and 
are invariant under this rotation, as may be verified by the character table, Table |. 

T' is generated by the rotations and C3 depicted in Fig. 1. Because of the double- 
valued nature of T' rotations, these elements actually have orders 4 and 6, respectively. For 
reasons to be described below, it turns out to be convenient to present explicit representa- 
tions (reps) for an element of order 4 (such as C 2 ) and one of order 3 (such as C 3 R). We 
label these elements as g 5 and g 9 , respectively;^] then T' is defined by the multiplication 
rules gl = g\ = 1, g % g\ = g^g 9 , and g 5 gg 1 g5 = #9#5#9- One may then show that each of the 
24 elements may be written uniquely in the canonical form ggglgl, where p = 0, ±1, and if 
q = or 2 then r = 0, while if q = ±1 then r = 0, ±1. 

The group T' is central to our model building since it is the smallest with 1-, 2-, and 
3-dimensional reps and the multiplication rule 2 © 2 = 3 © 1. T" models therefore allow for 
flavons that perform the same roles as a , S a b and A a b in the U(2) model. The only other 
24-element group that has reps of the same dimensions is the octahedral group O (which is 
isomorphic to S4). In this case, however, the product of two doublet reps does not contain 
a triplet, and the analogy to U(2) is lost. 

More specifically, T' has three singlets 1° and 1 ± , three doublets, 2° and 2 ± , and one 
triplet, 3. The triality superscript provides a concise way of stating the multiplication rules 
for these reps: With the identification of ± as ±1, trialities add under addition modulo 
three, and the following rules hold: 

1 © R = R ® 1 = R for any rep R, 2 © 2 = 3 © 1, , , 

2 ® 3 = 3 ® 2 = 2° © 2+ © 2-, 3 © 3 = 3 © 3 © 1° © 1+ © 1~. ^ } 



For a review of basic terms of discrete group theory, see Ref. p2[, Appendix A 



3 The element labels are chosen to coincide with those of Thomas and Wood, Ref. [23], where T 1 
is seen to be isomorphic to SL2(F3), the group of two-by-two unimodular matrices whose elements 
are added and multiplied as integers modulo 3. 
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Note that trialities flip sign under Hermitian conjugation, so that 2 + ® 2 = 3 © 1° while 
(2+)t®2" = 3 1+. 

The multiplication of X" representations may be made explicit by the use of Clebsch- 
Gordan matrices. For example, let the fields x an d ip be column vectors that transform as 
2 + and 2 under T', respectively. From the rules above, we know that the product of these 
reps contains a trivial singlet, the 1°, but it is not immediately clear how to construct this 
representation out of the given fields. Formally, we seek a matrix M such that the product 

X T Mij - X T M^ (3.2) 

under the transformations x ~* R{o)X an d V* ~~ * R{o)Xi where R is a two-dimensional 
matrix rep, and g runs over all elements of the group. From our earlier discussion, it is only 
necessary that we consider transformations associated with the defining elements, #5 and gg, 
to solve for the form of M; in the present case, one finds that M is proportional to the Pauli 
matrix cr 2 . This algebraic procedure is easily generalized to products of other reps. Explicit 
matrix representations for the generating elements g§ and gg, as well as the complete set 
of Clebsch-Gordan matrices for combining T' reps are provided in Appendix A. The reader 
should keep in mind that these Clebsch-Gordan matrices must be taken into account if one 
is to reproduce the Yukawa textures presented later in this paper. For example, without the 
factor of o"2, one might not realize that a vev in the first component of x couples only to the 
second component of ijj. 

As mentioned in the Introduction, we also require that our discrete flavor symmetry 
contain a subgroup that rotates first-generation matter fields by a phase. This subgroup 
plays the same role as the intermediate U(l) symmetry in the U(2) model, and must forbid all 
entries in the first row and column of each Yukawa matrix. The smallest discrete subgroup 
that one might consider is a Z 2 that flips the sign of all first generation matter fields. 
Unfortunately, such a transformation leaves the 11 entry of each Yukawa matrix invariant 
(two sign flips), so that the up and down quarks could, in principle, acquire masses that are 
too large. A Z3 phase rotation, on the other hand, does not lead to the same problem, and 
a Z 3 subgroup of T' is generated by the element gg defined previously. From Appendix A, 
we see that the two-dimensional representation matrices for the element gg are given by 

^°)=(; 2 J), *( 2+ )=(i£)> »( 2 >(;!)> ( 3 - 3 ) 

where 77 = exp(27ri/3). If matter fields of the first two generations are assigned to the 2~ 
rep, one then obtains the desired phase rotation under the Z 3 subgroup. This observation 
is at the heart of the global T" model presented in Section VII.] 



As we see below, however, models in which T' is free of discrete gauge anomalies are 
much easier to construct if matter fields are assigned to the 2° rep instead. In this case, 
let us consider extending the flavor symmetry group to T' x Z3. We identify a new triality 



4 One can also imagine models in which the symmetry group breaks to a non-Abelian subgroup; 
however, in this case the simple rephasing of multiplet components under the subgroup is not 
guaranteed. 
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index 0, + and — with the Z% phase rotations 1, n, and n 2 , respectively. Like the T' indices, 
the Z 3 trialities also combine via addition modulo 3. Reps of X" x Z 3 are denoted by affixing 
this additional triality as a superscript, e.g., 2 +_ . We now identify the desired intermediate 
symmetry as the diagonal subgroup of the original Z 3 , generated by the element g 9 , and 
the new Z% factor. We call this subgroup Z® henceforth. It is easy to see that the rep 2 0_ 
transforms under Z® by the matrix 

(3.4) 




which is simply the product of gg(2 ) and n . The matter field assignments 2 © 1 , and 
the breaking pattern T 1 x Z 3 — > Z® — > nothing are at the heart of the minimal flavor model 
discussed in the next section. It is worth pointing out that the reps l 00 , 1 H , 1 h , 2 0_ , 
2 ++ and 2~° are special in that these singlet reps and the second component of the doublets 
remain invariant under Z%. Thus any 2 © 1 combination of these reps is potentially useful 
in building models with U(2)-like textures. 

Finally, we return to the issue of anomaly cancellation. We pointed out in Section J I 
that consistency of a discrete gauge symmetry at low energies only requires the cancellation 
of anomalies that (1) involve the unbroken non-Abelian continuous gauge groups, and (2) 
are linear in a continuous group in which the discrete group is embedded. If we embed 
T' in SU(2), then these constraints are satisfied automatically, providing that the particle 
content of a given model fills complete SU(2) representations. Let us therefore consider the 
embedding of T' in SU(2) in more detail. 

The group SU(2) has one rep of each nonnegative integral dimension n (the spin (n— 1)/2 
rep), while T' has only singlet, doublet, and triplet reps. It must be the case that large 
SU(2) reps break up into a number of T' reps with the same total dimension. To see this 
decomposition, consider the characteristic polynomial of matrices in each of the T' reps for 
any two rotations that generate the full group. The same can be done for the full SU(2) 
group restricted to the particular Euler angles that give T'. Then a large rep matrix of 
SU(2) is block-diagonalizable into smaller blocks corresponding to rep matrices of T"; in 
particular, the characteristic polynomial of the SU(2) matrix is the product of those of the 
T' matrices. It is then possible to extract which T' reps appear in a given SU(2) rep, as 
well as their multiplicities. The results of this decomposition are summarized in Table ||. 
There we see that the 1°, 2°, and 3 reps of T' correspond to the complete 1, 2, and 3 
reps of SU(2). It follows, for example, that T' is non-anomalous in all models utilizing the 
2 0_ © l 00 representation structure for the matter fields (with Higgs fields as singlets). Note 
that there is no meaningful low-energy constraint on the Z% charges since Abelian factors 
may be embedded at high energies in U(l) gauge groups whose anomalies are cancelled by 
the Green-Schwarz mechanism fbl . 



IV. A MINIMAL MODEL 

In this section we review the minimal T" x Z 3 model presented in Ref. [12], which we 



study in quantitative detail in Section |V|. The three generations of matter fields are assigned 
to the T' x Z 3 reps 2°~ © l 00 while the Higgs fields Hjj,d are taken to be pure Gf singlets. 
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Given these assignments, it is easy to obtain the transformation properties of the Yukawa 
matrices, 



Y, 



U,D,L 



[3~ © I "] 


[2°+] 


[2 0+ ] 





(4.1) 



Equation ( j4.1|) indicates the flavon reps needed to construct the fermion mass matrices, 
namely, 1 0_ , 2 0+ , and 3~, which we call A, <f), and S, respectively. Once these flavons 
acquire vevs, the flavor group is broken. We are interested in a two-step breaking 

T' ® Z 3 nothing, (4.2) 

where e' < e again represent ratios of flavon vevs to the scale Mf. Since we have chosen a 
'special' doublet rep for the first two generations, which transforms as diag{i], 1} under Z%, 
only the 22, 23, and 32 entries of the Yukawa matrices may develop vevs of 0(e) originating 
from vevs in S and <fi. The symmetry is then broken by a 1 0_ vev of O(e'). The Clebsch- 
Gordan coefficient that couples a 1 0_ to two 2°~ doublets is proportional to <r 2 , so the e' 
appears in an antisymmetric matrix. These considerations yield the textures 



Yu,d,l ~ -e' e e , (4.3) 




where 0(1) coefficients have been omitted. Since the 1 0_ and 3~ flavon vevs appear as anti- 
symmetric and symmetric matrices, respectively, all features of the grand unified extension 
of the U(2) model are obtained here, assuming the same GUT transformation properties are 
assigned to (ft, S, and A. One can also show readily that the squark and slepton mass-squared 
matrices are the same as in the U(2) model. 

It was shown in Ref. [|T^] that this simple model can be extended to describe the observed 
deficit of solar and atmospheric neutrinos. Models for lepton masses were constructed both 
with and without the assumption of SU(5) unification. The latter possibility is of interest, 
for example, if one is only concerned with explaining flavor physics of the lepton sector, and 
is provided for completeness. In either case, the proposed extensions yield viable neutrino 
textures with naturally large mixing between the second and third generations. Moreover, 
these extensions do not alter the charged fermion textures, so that all the relations between 
masses and mixing angles in the U(2) model are also predictions of T' x Z%. We now review 
the two cases considered in Ref. |T2| . 

Case I: Here we do not assume grand unification, so that all flavons are SU(5) singlets. 
We introduce three generations of right-handed neutrinos transforming as 

v R ~ 2°- © 1-+ (4.4) 

Note that this representation choice differs from that of the other matter fields only in the 
third generation. Since Ur are singlets under the standard model gauge groups, introducing 
a 1 field by itself creates no anomaly problems. The neutrino Dirac and Majorana mass 
matrices then allow flavons that do not contribute to the charged fermion mass matrices. 
Their transformation properties are given by 
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M, 



LR 



[3- © I "] 


[2«] 


[2 0+ ] 





M, 



RR 



[31 


[2" 


hOj 


[2+0] 


[1- 





(4.5) 



Note that one obtains the same triplet and nontrivial singlet in the upper two-by-two block 
as in the charged fermion mass matrices, as well as one of the same flavon doublets, the 
2°+; the rep 1°~ is not present in Mrr, since Majorana mass matrices are symmetric. In 
addition we obtain the reps 2+°, l+~, and which did not appear in Eq. (|4.1| ). New 
flavon fields can now be introduced with these transformation properties, and their effects 
on the neutrino physics explored. Let us introduce a single^ new flavon $ v transforming as 
a 2+° and with a vev 




(4.6) 



where a 2 is the Clebsch that couples the two doublets to This new flavon is the only 
extension we make to the model in order to describe the neutrino phenomenology. After 
introducing (f) u , the neutrino Dirac and Majorana mass matrices read 



he 1 hr 2 e> 
M LR « | -l x e' l 2 e hne | (H v ) 
l A e 



M, 



RR 



/ /2 / 

/ r^r 2 e r 4 riee 



r A riee 



V r 2 e' 




(4.7) 



where is the right-handed neutrino mass scale, and we have parameterized the 0(1) 
coefficients. Furthermore, the charged lepton Yukawa matrix including 0(1) coefficients 
reads 



Yr. 



I cie' 
-c x e' 3c 2 e c 3 e 
V c 4 e c 5 



(4.8) 



The factor of 3 in the 22 entry is simply assumed at present, but originates from the Georgi- 
Jarlskog mechanism in the grand unified case considered next. 

The left-handed Majorana mass matrix Mn follows from the seesaw mechanism 



M, 



LL 



M LR M^ R M T LR 



(4.9) 



which yields 



M. 



LL 



( (e'/e) 2 e'/e e'/e 

e'/e 1 1 
V ^/e 1 1 



(Hu) 2 e 



A 



R 



(4.10) 



where we have suppressed the 0(1) coefficients. It is clear by inspection that we naturally 
obtain large mixing between second- and third-generation neutrinos. It is also important to 
point out that the two eigenvalues of Eq. (|4.10| ) that appear to be of 0(1) depend sensitively 



5 Assuming more than one (f) u leads to the same qualitative results. 
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on the products of a large number of order one coefficients. It is easy to obtain a hierarchy of 
order 10 in the two largest mass eigenvalues, without allowing any of the coefficients defined 
in Eqs. ( |4.7|) — (|4.8|) to deviate from unity by more than a factor of 2. This comment is 



important in understanding how the reasonable coefficient choices given in Ref. [12| account 
for the differing mass scales associated with atmospheric and solar neutrino oscillations. 

In order to determine neutrino oscillation parameters precisely one needs to compute the 
neutrino CKM matrix. If M LL and Y L are diagonalized by M LL = WMl L W\ Y L = UiXgUk, 
then 

V = U ] L W. (4.11) 



We parameterize this matrix as in Ref. 



V 



C12C13 C13S12 S13 



-C23Si2e l</> - C12S13S23 ci 2 C22,e l<t> - S12S13S23 C13S23 I ■ (4.12) 

V s 23 s 12 e — C12C23S13 — Ci2S23e l<?i — C23S12S13 C13C23 



where Cij(sij) stands for cos 0y( sin 6ij). Then one finds 

sm2 < 2 ^) = 4 (v|tW' (413) 

S in 2 (2fe 3 ) = 4 Jfi . (4.14) 

(V 2 3 + V33J 

The observed atmospheric neutrino fluxes may be explained by v^-v? mixing if sin 2 2# 2 3 ^ 
0.8 and 10~ 3 < Am^ < 10~ 2 , while the solar neutrino deficit may be accommodated by 
v e — mixing assuming the small-angle MSW solution 2 x 10~ 3 < sin 2 20 12 ^ 10 -2 for 
4 x 10 -6 % Am 2 2 ^ 10~ 5 , where all squared masses are given in eV 2 [^,^]. These regions 



of parameter space are the ones obtained most naturally from our models.^ Since is not 
determined from symmetry considerations, it is only necessary to reproduce Am| 3 /Am 2 2 . 
In Ref. [0] a choice for the 0(1) coefficients can be found that yields neutrino mass ratios 
and mixing angles falling within the desired ranges given above. 

Case II: Here we assume SU(5) unification and that the flavons transform nontrivially 
under the GUT group, namely, A ~ 1, S ~ 75, <p ~ 1, and S ~ 24. Note that since H ~ 5, 
the products SH and AH transform as a 45 and 5, respectively, ultimately providing a factor 
of 3 enhancement in the 22 entry of (the Georgi-Jarlskog mechanism). In addition, two 
2 +0 doublets are introduced, <$> v \ and 4> u2 , since the texture obtained for the neutrino masses 
by adding only one extra doublet is not viable. Both doublets (f) v have vevs of the form 
displayed in Eq. ( |4.6| ). As before, the presence of these two new doublets does not alter the 
form of any charged fermion Yukawa texture. 

The neutrino Dirac and Majorana mass matrices now take the forms 



6 The experimental ranges for neutrino mixing parameters follow from a two-neutrino mixing 
approximation which is valid only if the mixing angle #13 < 15° This condition is satisfied in 
all our models. 
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A 



(4.15) 



while the charged lepton mass matrix is the same as in Eq. ( |4.8| ) . Using Eq. (|4.9|) one obtains 
the texture: 



M, 



LL 



e '/ € 1 1 / Ab 



(4.16) 



Again, a viable set of 0(1) coefficients may be found in Ref. [|i~2 



While the texture in Eq. (|4.16j ) appears to be the same as the one in Eq. ([Oi]) (up to 
an overall factor of e), there is in fact an important difference: the 0(1) entries in Eq. ( [4.16| ) 
have a vanishing determinant at lowest order. The ratio of the two largest eigenvalues 
are therefore determined by higher order corrections, which must be taken into account to 
obtain the correct numerical results.[] While the zero determinant is lifted at 0(e) in the 
superpotential, it is interesting that, in this particular case, a comparable correction comes 
from D-terms that alter the canonical form of the neutrino kinetic energy 



d 4 e[ulu L + v\Bv L \ 



(4.17) 



Here B is a Hermitian matrix that depends on the flavons in the model. The kinetic terms 
may be put back into canonical form by the superfield redefinition 
(1 — B /2)z/£. This in turn leads to a correction to Mil, 



M LL ^ M LL -1/2{B,M LL }. 
Numerically, it is only necessary that we retain the largest elements of B 

(■ ■ ■ 



^n^Bv L « 

(4.18) 



B 



ae 



(4.19) 



\ • ae 



which also leads to an 0(e) correction to the determinant discussed above. The parameter 
a is included in the quantitative analysis of the model presented in the next section. 



V. NUMERICAL ANALYSIS 



The numerical check of the unified T' x Z3 model presented in [0 relied on two assump- 
tions. The first is that there exist 0(1) coefficients q, di, and Ui for the charged fermion 
Yukawa matrices that, when combined with the particular choice of neutrino Yukawa pa- 
rameters U and rj, yield charged fermion mass eigenvalues and mixing angles in agreement 



7 In fact, the analysis made for the model in Case I included higher order terms, which did not 
contribute in any significant way. 
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with the values observed. This should not be a problem since the textures of the X" x Z; 



3 



model for the charged fermions agree completely with those of the U(2) model [10|, in which 



all of these observables are accommodated in detailed fits. Second, the textures as written 
in the last section are defined at the scale M GUT ~ 2 x 10 16 GeV, while the observables 
are of course measured below the electroweak scale. A truly meaningful fit requires running 
the gauge and Yukawa couplings over this range. While the textures renormalized at Mqut 
and m t should not differ wildly in form, a global fit is required to properly compare the 
predictions of our model to the experimental data. The purpose of this section is to report 
on the necessary steps in these fits and the numerical results. 

In order to study the renormalization of gauge and Yukawa couplings, we run the one- 
loop renormalization group equations (RGE's) of the MSSM p7| from Mgut down to the 
electroweak scale taken to be m t = 175 GeV. This analysis does not include two-loop 
corrections nor threshold effects at either end of the spectrum. In particular, this approach 
does not differentiate between the scales Mf, eMf M GU t, « e'M GUT /e, and A R « 
eMcuT-Q In any case, both the two-loop and threshold effects are formally of subleading 
order, and therefore are taken into account by permitting theoretical uncertainties in the 
gauge and Yukawa couplings of 0(l/167r 2 ) ~ 1%. 

Values of the gauge couplings at M GU t are obtained by starting with the precision values 
extracted at the scale M z P^fl , 

a^(M z ) = 59.99 ±0.04, 
oq\M z ) = 29.57 ±0.03, 

a^{M z ) = 8.40 ±0.13. (5.1) 

The gauge couplings are run from M z to m t using the one-loop Standard Model (SM) RGE's, 
and then from m t to M GU t using the one-loop MSSM RGE's.[] The GUT scale couplings 
are taken directly from the textures of Eqs. (|1.6|) , fll.9Q , (|4.8| ), and (|4.15| ), given numerical 
values for the dimensionless coefficients q, di, k, r^, Ui, and a (collectively fcj), and for e, e', 
p, and £. The Yukawa matrices are then run down to m t and diagonalized.0 



8 Notice that Ar ~ eM G uT yields the appropriate mass scale in Eq. ( [4.16 ) for atmospheric neutrino 
oscillations. 

9 It should be pointed out that, while the SM RGE's make use of the MS scheme, the MSSM 
RGE's in Ref. [27| make use of the DR scheme [^], which differ at the matching scale (mt by our 
choice) by an amount 

4tt 47T 1 

where C A = {0, 2, 3} for i = 1, 2, 3. 

10 The RGE's are integrated by means of the Runge-Kutta method with adaptive stepsize control 
[ ]30[ j . The results of this method were cross-checked against the results of using Richardson ex- 
trapolation with Bulirsch-Stoer stepping [3C] and were found to agree to the limits of the expected 
accuracy of either solution. 
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Experimental values for the low-energy Yukawa couplings are extracted from the physical 
masses and mixing angles compiled by the Particle Data Group |28| , where entries of Yjj are 
obtained by dividing quark masses by v sin/^/v^ and those of Y D L by dividing quark and 
lepton masses by v cos P/y/2 , where v = 246 GeV. 

The experimental uncertainties on the observables (or estimates for the quark masses) 
used in the fits are either those appearing in Ref. or 1% of the central value, whichever is 
larger; since the lepton masses are measured with extraordinary precision, they are sensitive 
to the two-loop RGE and threshold corrections that we have ignored. 

The RGE for the neutrino Majorana mass matrix Mll was computed in Ref. BTJ and is 



included here in order to complete the RGE evolution for all observables. The low-energy 
neutrino observables are taken to be 



Am\ 

'-12 



100 < — < 2500, 
Ami 



sin 2 26 23 > 0.8, 

2 x 10~ 3 < sin 2 20i2 < 10~ 2 . (5.3) 

For the sake of having meaningful uncertainties, a parameter whose lower bound is much 
smaller than its upper bound is converted into its logarithm. Instead of Eq. ( |5.3| ), we use 

sin 2 2# 23 = 0.9 ±0.1, 
In (sin 2 29 12 ) = -5.41 ± 0.80. (5.4) 

Summarizing to this point, we have discussed the details of how inputs consisting of the 
gauge couplings at Mz and Yukawa matrix parameters at a high scale are manipulated using 
one-loop RGE's to produce output values for fermion masses and mixing angles observed at 
low energy. Of course, the salient question is whether one can find a choice of parameters ki, 
where all of these coefficients are 0(1), and yet the output quantities are all in agreement 
with their observed values.^] This is accomplished through a x 2 minimization; thus, the 
complete simulation consists of choosing a set of parameters ki (relevant at Mqut), running 
the RGE's down to m t , and comparing with observation to compute a figure of merit, 
X 2 . If X 2 is too large, the parameters ki are adjusted and the procedure is repeated until 
convergence of x 2 to a minimum is achieved. 

The x 2 function assumes a somewhat nonstandard form. Fermion masses and mixing 
angles are converted to Yukawa couplings k^ xpt ± Aki, and contribute an amount 

/ , cxpt _ , \ 2 

A * 2 = [- L M A ) < 5 ' 5 > 



11 We also allow for variation of the parameters e, e', p, and £ by hand, but do not minimize 
with respect to them. Changes in these parameters are equivalent to redefinitions of the 0(1) 
coefficients, so that they merely set the scale for the other parameters of the fit. 
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to x 2 ■, as usual. There are 15 observables (6 quark masses, 3 quark CKM elements [since 
CP violation is neglected], 3 lepton masses, 2 neutrino mixing angles, and 1 neutrino mass 
ratio) and 26 parameters kf, on the surface, it seems that the fit is always under-constrained. 
However, our demand that the parameters ki lie near unity imposes additional restrictions, 
which we include by adding terms to x 2 °f t ne form 

A^(l^) 2 (5,, 

for each i. Thus, the parameters ki are effectively no longer free, but are to be treated 
analogously to pieces of data, each of which contributes one unit to x 2 if it is as large as 
3 or as small as 1/3. The particular choice of 3 for this purpose is, of course, a matter 
of taste. In effect, the inclusion of such terms renders the parameters ki no longer as true 
degrees of freedom. On the other hand, they are not true pieces of data either, since a value 
of say, ki = 0.8 is just as valid as a value of — 1.1 for our purposes. Thus, the value of 
Xmin determining a 'good' fit is 15, since there are 15 pieces of true data and effectively no 
unconstrained fit parameters. 

The numerical minimization is carried out using the MINUIT minimization package. 
As a cross check, minimization using Powell's direction set method [fj0| is carried out to 



make sure that the same minimum is achieved. Since the topography of the x 2 function is 
complicated due to the numerous parameters involved, it is important to try a number of 
initial choices for the input parameters ki in order to have confidence that the minimum 
obtained is close to global. Once convergence is achieved, a parabolic minimum is assumed 
and a Hessian matrix is computed in order to gauge uncertainties of the parameters. 

Detailed numerical fits show that it is not difficult to find parameters ki that satisfy the 
constraint Xmm < 15. However, in the T'xZ 3 model, the ratio mb/m t must be accommodated 
either by a small value of £ or a large value of tan j3. For definiteness, we choose tan j3 = 2 as 
a representative value, and find a best fit with Xmm °f 2.77. The complete set of parameters 
is listed in Table [TTT| and a comparison to data appears in Table [TV]. Note especially that 
the parameters e, e', and p are somewhat larger (a factor of 2 or more) than their values in 
the U(2) model of Ref. |9]], where neutrino physics was not considered. From the excellent 
X 2 , one concludes that the T' x Z 3 model has little difficulty satisfying all of the required 
constraints including the naturalness of the coefficients, allowing for the small parameter £ 
that distinguishes the scale of Y v from Y D L . 

While we have seen that the minimal scenario is extremely successful at reproducing 
fermion masses and mixing angles, there are nonetheless a number of interesting variant 
models based on T' symmetry. We explore these models in the next three sections. 



VI. SU(5)xU(2) WITH NEITHER SU(5) NOR U(2) 

As discussed in the Introduction, the U(2) model must be embedded in a grand unified 
theory to reproduce all of the observed quark mass hierarchies. In this section we present a 
model that does exactly the same, without the need for a GUT, by extending the discrete 
gauged flavor group to T' x Z 6 . We show that this model explains the ratio mb/mt, which 
is merely parameterized in the U(2) model (and in our other T' models). Before presenting 
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the model we comment on notation. As before, we use the triality superscripts +, — , and 
for the different representations of T' . For the Z§ reps we now introduce the indices 
% = 0, 1, . . . , 5, which combine through addition modulo 6. For example, 2 +4 © 1 +2 = 2 _0 , 
etc. Since Z§ is isomorphic to Z 3 x Z 2 , one may view the new flavor symmetry as a Z 2 
extension of the T' x Z% flavor group defined in the model of Section [IV]; denoting the Z 2 
reps by + and — , one identifies 





+ 



+ 




1 

2 
3 
4 

5 



That is, the Z 6 charge is 2 x (Z 3 charge) + 3 x [Z 2 charge) modulo 6. In the remainder of 
this section we use the more compact T' x Z 6 notation. 
The three generations of matter fields transform as 



Q,U,D~2 



04 



I 00 



(6.1) 



L ~ 2 04 © 1 



E ~ 2 + © 1 



(6.2) 
(6.3) 



,04 



l+l 



(6.4) 



The matter fields have transformation properties that differ from those in our previous 
models, and in particular, the electroweak doublet leptons are no longer anomaly free by 
themselves. The third-generation L field is assigned to a nontrivial T' singlet, the 1 + , which 
does not form a complete SU(2) representation. Given the discussion in Section |I|, the 
T' SU(2)^ anomaly is not automatically cancelled. However, we remedy this problem by 
assigning non-trivial transformation properties to the Higgs fields: 



u 



l 00 



H 



D 



(6.5) 



The fields Hp and L 3 are both electroweak doublets and, as far as the non-Abelian quantum 
numbers are concerned, form a vector-like pair when Hd is a 1~ under T' . The remaining 
fields, E and vr, do not transform under any unbroken non-Abelian continuous gauge groups 
and thus their T' x Z$ quantum numbers may be assigned freely. 

In order to break the flavor symmetry and obtain the fermion mass matrices we introduce 
the following flavons: 



5-3°, A 
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(6.6) 
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(6.7) 



In addition to these flavon fields, we introduce two more in the neutrino sector of the theory. 
Their transformation properties are such that they do not alter the form of the charged 
fermion Yukawa textures: 



+1 



(6.8) 



Together with Ur, these fields are the only ones that transform nontrivially under the Z 2 
subgroup of Z 6 (i.e., the only ones with odd Z 6 charges). Again, we are interested in a 
two-step breaking: 



T' x Z. 
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(6.9) 



where Z® is precisely the same subgroup as in the minimal T' x Z% model. Thus, by the 
same arguments presented in Section |IV|, we obtain the following patterns of vevs: 
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(6.10) 
(6.11) 
(6.12) 



Unlike the minimal model described in the previous two sections, the flavons here contribute 
to the Yukawa matrices in some cases only at quadratic order 
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(6.15) 



We see that the flavons A and A' appear in precisely the right way to recover approximate 
SU(5) x U(2) textures for Yd and Yl, with an additional overall factor of e. The only 
difference is a relatively uninteresting e' entry in the 13 and 31 elements of Yl- Notice that 
the vev of the E field has been replaced by (A) in Eq. ( 6.13| ). Thus, all important features 
of the SU(5) x U(2) model are reproduced. 

Note that the ratio mb/m t , which is experimentally observed to be in the range 0.023 < 
THi/rrit ^ 0.026, is predicted to be of order e ~ 0.02 for tan/3 ~ 0(1), as can be seen from 
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the ratio of the 33 entries in Yjj and Yd- This is promising since tan/5 ~ 0(1) is the naive 
expectation if the Higgs potential is not fine-tuned. 

Before proceeding to the analysis of the neutrino sector, a few comments are warranted 
on the possible supersymmetric contributions to FCNC's in this model. As mentioned in 
the Introduction, scalar superpartners of the first two generations are exactly degenerate in 
our models when the flavor symmetry is unbroken. The amount of scalar nondegeneracy 
at low energies is determined by the order at which flavons contribute to the scalar mass 
matrices. In the minimal model, the flavons contribute quadratically to the scalar masses 
of the first two generations, as a consequence of the flavons' nontrivial Z 3 charges. The 
scalar mass-squared matrices of the U(2) model are then reproduced. In the current model, 
however, the flavon S may contribute linearly, since 3° is in the product of (2 04 )^ <g> (2 04 ). 
The important point is that this effect provides an 0(e) correction to the diagonal entries 
of the scalar mass matrices. In the fermion mass-eigenstate basis, a Cabibbo-like rotation 
O c ~ e'/e leads to 12 entries in the scalar mass matrices of order e'rh^, where is an 
average scalar mass, and e' ~ 0.004. Taking into account uncertainty in 0(1) coefficients, 
this result is in marginal agreement with the bounds from CP-conserving flavor-changing 



processes, assuming superpartner masses less than a TeV While bounds from CP- 

violating precesses are generically stronger, the 0(1) coefficients have unknown phases that 
one may simply choose in order to avoid these bounds. Without a firm understanding of the 
origin of CP violation, saying more about these phases entails a degree of speculation that 
we choose to avoid. Of course, if scalar superpartners are heavy (as in the 'more minimal 
MSSM' ||33]| ) or flavor universal (as in gauge mediation anomaly mediation [@,|J, or 
Scherk-Swartz mechanism M) the current T' model is completely safe. 

Next we examine the neutrino sector of the model. Given the transformation properties 
of Ur, we calculate the neutrino Dirac and Majorana mass matrices 



M, 



LR 



[3 4 © l 04 ] 




[2- 4 ] 





AS 4 


- AA -f 


-0 2 


A(j> v 


A' 4 


' + 


K 


A v 



(H; 



V 



he' hne' 
-he' he l 2 r 3 e \ e(H v } , (0.10) 
he he h 



M. 



RR 



[3 4 ] 




[2- 1 ] 


[1 +4 ] , 



AS 


A<j) u 


A<f> v 


A 



A 



R 








ne' 

r 2 e r 3 e 



eA 



R 



(6.17) 



where and h are 0(1) coefficients. To leading order, the seesaw mechanism gives 
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(6.18) 



Note that the texture in Eq. ( |6.18|) is not changed if higher-order corrections are included 
that lift the zeroes in Eqs. flb.lfip — ( |6.17| ). Following the same procedure as before, we 
diagonalize Mll and Yl and extract the neutrino masses and mixing angles. A global fit 
of the parameters in this model can in principle be done; however we just present a viable 
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set of parameters for simplicity. Using the set of values for the 0(1) coefficients in Mn ( 
r 1 ,...,r 4 ,/ 1 ,...,/ 6 ) = (1.0, 1.0, 1.0, -1.0, 1.2, 1.2, 1.3, -1.0, -2.0, 1.0) and assuming all 
coefficients in Yl are 1.0 except that of the 22 entry, which we set to 3.0, we obtain 



Am|3 
Am 2 12 



125, 



sin 2 26 



12 



3.5 x 10" 



sin 2 26 
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0.88. 



(6.19) 



This agrees with the allowed ranges described in the previous sections. It is worth mentioning 
that the texture Eq. (|6.18|) is the same as obtained in Ref . [ 34j] , and thus the claim in Ref . |35 
that this texture cannot account for solar neutrino oscillations is not correct. 



VII. A GLOBAL T' MODEL 

As pointed out in the Introduction, it is not possible to construct a realistic super- 
symmetric model with a continuous SU(2) flavor symmetry if scalar universality is not 
assumed. The argument is straightforward: The left- and right-handed up quark fields must 
be embedded in 2©1 representations to maintain the heaviness of the top quark, as well 
as degeneracy of squarks of the first two generations. Given this assignment, the coupling 
Q a U b e a bH u is allowed by the unbroken flavor symmetry, which implies the unacceptable re- 
lation m u = m c « m t . The T' model below demonstrates that discrete subgroups of SU(2) 
are viable for building models of fermion masses, although they are more dangerous than 
models with additional Abelian factors, as far as supersymmetric FCNC processes are con- 
cerned. We first present the model, and then explain how it evades the problem described 
above. 

The crucial feature that allows one to build a successful T" x model is the existence 
of a doublet representation 2 0_ , whose first generation component alone rotates by a phase 
under the Z® subgroup. This choice is unique in models where T' is a discrete gauge 
symmetry, since the 2° rep is the only doublet that fills a complete SU(2) representation if 
we embed T' in SU(2). The 4 of SU(2) decomposes into the reps 2 + and 2~, which implies 
that each is separately anomalous. While it might still be possible to construct models 
involving anomaly-free combinations of 2 + and 2~ reps, we have found no examples that are 
particularly compelling. On the other hand, if T' is assumed to be a global symmetry, then 
matter fields can be assigned to any of the doublet representations freely. This provides an 
opportunity for constructing economical models, as we now demonstrate. 

Consider the Z 3 subgroup of T' generated by the element gg that acts on the 2° rep as 
diag{r] 2 ,r]}, with rj defined as in Section |Hl[ In the 2~ rep, this element takes the form 
diag{r], 1}, which we identify as the desired phase rotation matrix for matter fields of the 
first two generations. Given our freedom to assign matter fields to any of the doublet reps 
in a global T' model, it is no longer necessary to extend the flavor symmetry by an Abelian 
factor in order to find a subgroup that forbids the order e 1 Yukawa entries. Thus, one is 
naturally led to the charge assignment 

V>~2-©1° for V = Q, U, D, L and E, (7.1) 

and Hjj,d ~ 1°, which yields 
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Y, 



U,D,L 



[3 



[2- 



(7.2) 



Introducing flavons, A, <ft and S transforming as 1 , 2 + , and 3, respectively, one reproduces 
the canonical U(2) textures assuming the breaking pattern 



r 



nothing, 



(7.3) 



together with the dynamical assumption that only the 1~ rep participates in the last step 
of symmetry breaking. The resulting textures are identical to those in our original model 
of Section [TV]. One difference, however, is that the S flavon in this model contributes to the 
squark mass matrices at first order in e, just as in the T' x Z$ model. However, this is not 



a concern for the same reasons discussed at length in Section [VI 



Turning to neutrino physics, recall that successful results were obtained in the T' x Z% 
model by altering the charge assignment of the third-generation right-handed neutrino field. 
Thus, we are motivated here to consider 



which implies 



M: 



LR 



M, 



RR 



[3] 



(7.4) 



(7.5) 



We identify the flavon <fi v with the representation 2~, which does not appear in any of the 
charged fermion Yukawa textures. However, there is an important difference between this 
model and the one discussed in Section [TV]: The third generation v R field transforms by a 
phase under the Z% subgroup, so that, for example, the 13 and 31 entries of M RR are left 
invariant under this intermediate symmetry. This implies an inversion in the hierarchy of 
vevs in the third row and column of Mrr. In the non-unified version of the model, it is 
somewhat remarkable that we still obtain a viable form for Mil'- 



he' hr^e 
M LR « | -l x e' l 2 e l 3 r 2 e' | (H v ) 
he 



M, 



RR 



/ ne 

r 3 e r 2 e' \ A R 
\ rie r 2 e' r 4 e' 



(7.6) 



e'/e 1 1 



Mll - | e'/e 1 1 | '-fLl, (7.7) 

Unfortunately, this result does not persist in the simplest unified version of the model, which 
includes additional suppression factors in the 22 entries of M LR and M RR . Fortunately, a 
simple modification of the flavon charge assignments in the unified theory allows us to 
recover the previous result. We introduce two (f) u flavons that transform differently under 
T'x SU(5): 

V ~(2- 24) , #,~(2- 1) . (7.8) 
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Furthermore, we assume the pattern of vevs 



(0.) ~(e)' ~ ( ) • (7 ' 9) 

This is consistent with the breaking pattern in Eq. (|7.3| ), but includes a dynamical as- 
sumption that the doublet <p' u does not participate in the first stage of sequential symmetry 
breaking and its second component acquires no vev.0 Since <p u transforms as an SU(5) 
adjoint, it can contribute directly to Mlr, but only to M RR if, for example, the adjoint 
flavon E is also present; the corresponding entries of Mrr are therefore suppressed by an 
additional factor of e: 

he' hr x e \ / r x e 2 

Mu<*\ -he' ke 2 hr 2 e' \ (H v ) , M RR « r 3 e 2 r 2 e' | A R . (7.10) 
/ 4 e / \ rit 2 r 2 e' r 4 e' 

The seesaw mechanism then yields 



/ (e'/e) 2 e'/e e'/e 

M LL ~ 

V ^/e 1 1 



(Hu) 2 e 



r'/r I I | (7.11) 



where we used the numerical fact that e' 2 /e 3 ~ 0(1). It is important to note that we 
have only displayed the contributions to Eq. ( [7.1 0| ) linear in 0, S and A, for convenience; 
quadratic and higher order corrections lift the zero entries of these textures, but do not 
change the result in Eq. ( |7.11p qualitatively. Note that Eq. ( [7.11] ) is the same successful 
texture obtained in our original T' x Z 3 model. 

Finally, we return to the no-go theorem presented at the beginning of this section. It 
is not possible to construct a realistic model with a continuous SU(2) flavor symmetry and 
2 © 1 rep structure because an unwanted flavor-invariant operator may be formed from the 
product of two doublet matter fields. In our global T' model we have the freedom to assign 
matter fields to new doublet representations whose products contain no trivial singlets, thus 
avoiding the problem. 



VIII. T WITH STERILE NEUTRINOS 

In this section we comment briefly on the possibility of four light neutrino species. Rather 
than investigating the (vast) space of possible models, we simply show how the results of a 
successful extension of the U(2) model with a sterile neutrino proposed by Hall and Weiner 
(HW) |35| can be reproduced with T 1 symmetry instead. 

Consider a U(2) model with all matter fields, including three generations of right-handed 
neutrinos, in 2 © 1 representations. Given the canonical pattern of flavon vevs, one obtains 
a right-handed neutrino mass matrix of the form 



We consistently assume that a flavon that transforms nontrivially under a subgroup Hi either 
acquires a vev of order the scale at which Hi is spontaneously broken, or acquires no vev at all. 
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M, 



RR 




(8.1) 



Since Mrr is symmetric, there is no contribution from the flavon A, leading to a singular 
matrix. It is important to emphasize that the zero entries of Eq. ( |S.1| ) are not lifted at 
any order in e and e' as a consequence of the holomorphicity of the superpotential. From 
consideration of the U(2) indices of the flavon fields (or alternatively their charges under 
a U(l) subgroup of U(2)), it is possible to show that any contribution to the vanishing 



entries of Eq. (|8.1| ) requires the complex conjugation of a flavon field, which is not allowed 
by unbroken supersymmetry. If the pattern of flavon vevs is not altered, the first-generation 
right-handed neutrino remains in the low-energy theory as a sterile neutrino. 

This sterile neutrino mixes with the second-generation left-handed neutrino at order e' 
in Mlr. After integrating out the two heavy right-handed neutrino flavors, one obtains a 
four-by-four neutrino mass matrix of the form 



(8.2) 





/ 


\ 


= 




ce'(Hu) 





[o ce'(Hu) 





where the three-by-three block has entries of order (Hjj) 2 /A^, which can be found in 
Ref. |35|. HW observe that the 24 and 42 entries of M< 4 ' are much larger than all others, 
leading naturally to maximal mixing between and the sterile neutrino. As it stands, 
however, both would have masses of order of the electroweak scale unless c is taken to be of 
O(10~ 8 ). To obtain a viable model, HW extend the flavor symmetry by an additional U(l) 
factor, under which all the right-handed neutrinos have charge +1. A charge —1 flavon is 
introduced with the vev ~ 10~ 8 , which breaks this symmetry weakly. One then finds 
that c ~ e^, while remains unchanged. 

The main obstacle to implementing this solution in a T" x Z 3 model with all matter fields 
assigned to 2°~ © l 00 reps is that higher-order corrections to the first row and column of 
Eq. ( |S.1| ) are not forbidden by holomorphicity; the complex conjugate of any non-trivial Z 3 
phase rotation is the same as its square. Thus, we are led to promote our Z 3 symmetry to 
a continuous U(1).0 The appropriate embedding is given by 



, 00 



-2 ■ 



^ 2 0+ — ► 2° x , S ~3 — ► 3 2 , A ~ 1°- — 
where the subscript indicates the U(l) charge. Assuming the breaking pattern 



T' x 17(1) 



7D e' 



nothing, 



(8.3) 



U) 



13 We could also promote Z3 to a much larger Z n that adequately suppresses corrections to the 
zero entries in Eq. ( |3.lD ; we leave this possibility implicit in our discussion. 
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we reproduce the textures of the U(2) model, including Eq. (|8.1|) , identically. The HW 
predictions for solar, atmospheric and LSND |36| neutrino oscillations are then recovered 
by extending the symmetry by an additional U(l) factor, implemented precisely as before. 
We are thus able to reproduce the results of Ref. |35| with the flavor symmetry T'xU(l) 2 . 
Although we find this model less compelling than the other three already discussed, it may 
be of some relevance if the LSND oscillation result is independently confirmed. 

IX. CONCLUSIONS 

We have shown in this paper how to reproduce the quark and charged lepton Yukawa 
textures of the U(2) model using a minimal non-Abelian discrete symmetry, the double 
tetrahedral group X'. The first model we discuss, based on the discrete gauge symmetry 
X' x Z3, not only successfully accommodates the observed charged fermion masses and CKM 
angles, but also accounts for solar (small-angle MSW) and atmospheric neutrino oscillations. 
In particular, a large u fi -v T mixing angle is predicted in the model, even though all charged 
fermion Yukawa textures are hierarchical. A global fit including neutrino parameters was 
performed in a grand unified version of the model, and results with extremely good \ 2 were 
obtained. 

In addition, two variant X" models were discussed. In the first, the flavor group was 
extended to X" x Z 6 , and all important features of the SU(5)xU(2) model were reproduced 
without the need for a field-theoretic unification. This model provided a successful prediction 
(with order-one uncertainty) of the bottom to top Yukawa coupling ratio, which is merely 
parameterized in the U(2) model and in the other X' models we discuss. The second variant 
theory was based on a global X' symmetry and demonstrates that the successful U(2) textures 
can be obtained without including an Abelian factor in the flavor group. In both variant 
models, large y^-v T mixing is predicted, and solutions to the solar and atmospheric neutrino 
problems are naturally obtained. 

It is worth pointing out that the viable neutrino textures predicted by our models are 
achieved without altering the predictive textures of the charged fermions, and without in- 
troducing sterile neutrinos. Interestingly, the solutions we present have no simple analogy 
in the U(2) model: the right-handed neutrino fields in our models do not fill complete U(2) 
representations. In particular, the third generation i/r transforms as a 1~, which forms only 
part of a 5 in U(2). Aside from the possibility of very nonminimal U(2) models (e.g. with 
seven generations of right-handed neutrinos), the desired neutrino X' reps do not naturally 
occur. The key advantage of discrete groups is that the large, phenomenologically unused 
representations of the continuous embedding group break up into sets of small phenomeno- 
logically useful representations of the discrete group. If discrete gauge symmetries arise as 
fundamental symmetries of nature, then we see from the example of T' that their richer 
representation structure makes it possible to construct simple and elegant models of flavor. 
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APPENDIX A: EXPLICIT DETAILS OF T 

As described in the text, the group T' is generated by the elements labeled g$ and g$. 
We begin by exhibiting explicit matrices representing these elements in each of the seven 
reps listed in Table I. The singlets are g 5 (l 0,± ) = 1, #9(1°) = 1, g9(l + ) = f], g 9 (l~) = r/ 2 , 
where r\ = exp(27ri/3). The doublets are 

05(2°*) = M l5 g 9 (2°)=r)M 2 , g 9 (2+) = r, 2 M 2 , g 9 (2~) = M 2 , (Al) 

where 



J_ / +i + v / 2e^/ 12 \ t , ( 7] 
and the triplet rep is generated by 
<fc(3) = I 



(A3) 



277 


2 V 2 \ 




(10 fl\ 


-1 




, <fo(3) = 


0?) 


2 V 2 


2 V 




v r? 2 / 



The Clebsch-Gordan (CG) coefficient matrices coupling an n^-plet x and an n y -plet y to 
form an n z -plet z consist of n z matrices of dimensions n x x n y satisfying the condition 

RlO i R y = Y,( R *)ii°v i = h---,n z , (A4) 
i=i 

where Ri denotes the group rotation R in rep i. In a perhaps more familiar notation, the 
CGs above may be written 



(Pi), 



f x y 



z 



(A5) 



Note from Eq. (|A5|) that the CG matrices for Ri ® R2 are simply the transposes of those 
for R 2 <8> Ri, and thus are omitted below. The coefficients c below indicate multiplicative 
constants arbitrary in the definition Eq. (|A4j) . The CG coefficients for two singlet reps or 
any rep with 1° are all unity; the remaining CGs for products involving singlets are 

l*i g, 2 h = 2 tl+h , with Ox = c(l 0), 2 = c(0 1). (A6) 
l+®3 = 3, with Ox = c(0 1), 2 = c(l 0), C 3 = c(0 1 0). (A7) 
1-03 = 3, with Ci = c(0 1 0), C 2 = c(0 1), 3 = c(l 0). (A8) 

Next, let 

M 3 = la-i)(°l), m, fi0 ^ 



2 V 7 v 1 J ' 4 \o / ' 
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Then 



2° <g> 2° D 3, 2 ± <g> 2 T D 3 : 

O l = cM 3 , 2 = cM 4 , 3 = cM 5 . 

2° <g> 2° D 1°, 2 ± <g> 2 T D 1° : 
O = cM 6 . 

2°®2 + D3, 2^2^33: 

O x = cM 5 , 2 = cM 3) O s = cM 4 . 

2° <g> 2 + D 1 + , 2" <g> 2~ D 1+ : 
O = cM 6 . 

2° <g> 2~ D 3, 2+ (g) 2 + D 3 : 

d = cM 4 , C 2 = cM 5 , C 3 = cM 3 . 

2° <g> 2~ D 1", 2+02+DT: 
C = cM 6 . 



The remaining combinations are: 

2 0>± <g> 3 D 2 0>± : 



_ / 1 \ / 1 - i 

Ol = c (o l + <0 ' ° 2 = c (-l 



2° <g> 3 D 2 + , 2 + <g> 3 D 2", 2^3 3 2°: 

/Ol \ 

° 1 = C (00 1 + iJ' ° 2 = c { -1 

2° <g> 3 D 2", 2+ (g> 3 D 2°, 2~ <g> 3 D 2+ : 
n ( o 1 \ / 1 - i 

Ol = c li + iooj' ° 2 = c (o -1 



3 ® 3 D 3 S ® 3 fl : 
/2 

Oi = ci 



\ 



-1 +C2 -1 

\o -l o / \o i o ; 
o -1 o\ /o -1 o 

2 = ci | -1 + c 2 1 | , 
2 / \ 
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(A19) 



3 <g> 3 D r 



o 



3®3Dl + : = c 



3(g) 3 D 1" : = 




(A20) 
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TABLE I. Character table of the double tetrahedral group T'. The phase 77 is exp(27T7/3). 



SU(2) rep multiplicity 


T' rep decomposition 


UN 


2iv{2°e2+e2-} 


12iV + l 


i°eiv{i ei + ©i-e3-3} 


12iV + 2 


2°0 2iV{2 o ©2+©2-} 


12iV + 3 


3©iV{l ©l + ©l-©3-3} 


12iV + 4 


{2+ © 2-} © 2N {2° © 2+ © 2-} 


12iV + 5 


{1+ © I' © 3} © N {1° © 1+ © 1" © 3 • 3} 


12iV + 6 


(2iV + l){2 ffi2 + ffi2-} 


12iV + 7 


{1° © 2 • 3} © N {1° © 1+ © 1" © 3 • 3} 


12iV + 8 


2°ffi(2iV + l){2 ffi2 + ffi2-} 


12iV + 9 


{1° © 1+ © 1" © 2 • 3} © N {1° © 1+ © 1" © 3 • 3} 


12iV+ 10 


{2+ © 2-} © (2N + 1) {2° © 2+ © 2~} 
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TABLE II. Decomposition of SU(2) reps into reps of T' . N is any nonnegative integer. 









e = 0.04, 


P = 


0.08, e' = 0.004 


, e = 


= 0.017 






Cl 


= -0.93 ±0.01 


di 


= +1.33 + 0.45 


h 


= +0.85 + 0.62 


n 


= +0.94 ± 0.84 


Ul 


= +0.92 + 0.31 


C2 


= -0.46 ±0.03 


d 2 


= -0.81 ±0.26 


h 


= -1.01 ± 1.11 




= +1.06 + 0.95 


U2 


= +1.48 ± 0.70 


C3 


= -1.02 ± 1.13 


d 3 


= +1.55 + 0.67 


h 


= -0.97 ±0.75 


<"3 


= +1.03 + 1.12 


U3 


= -0.90 ±0.91 


C4 


= -1.03 ±1.15 




= +1.14 + 1.33 


k 


= -1.09 ± 1.04 


r 4 


= -1.07 + 1.05 




= +1.07 + 1.21 


C5 


= -0.90 ±0.01 


d 5 


= -1.29 + 0.12 


h 


= -1.11 ±0.79 


<"5 


= -0.97 + 1.03 


u 5 


= +1.84 + 0.95 


a - 


= +0.98 ± 1.06 



















TABLE III. Best fit parameters for the T' x Z 3 model with tan/3 = 2. The minimum x 2 = 2.77. 



31 



Observable 


Expt. value 


Fit value 


m u 


(3.3 ±1.8) x 1(T 3 


3.5 x 10~ 3 


m d 


(6.0 ± 3.0) x 10~ 3 


4.0 x 10~ 3 


m s 


0.155 ±0.055 


0.136 


m c 


1.25 ±0.15 


1.24 


m b 


4.25 ±0.15 


4.25 


m t 


173.8 ±5.2 


170.4 


m e 


(5.11 ± 1%) x 10- 4 


5.11 x 10~ 4 




0.106 ± 1% 


0.106 


m T 


1.78 ± 1% 


1.78 


\Vus\ 


0.221 ±0.004 


0.221 


\Vub\ 


(3.1 ± 1.4) x 10~ 3 


2.3 x 10~ 3 


\Vcb\ 


(3.9 ±0.3) x 10~ 2 


3.9 x 10~ 2 


Amf 3 /Amf 2 


100 - 2500 


526 


In (Am| 3 /Amy 
sin 2 26»i2 


6.22 ± 1.61 


6.27 


2 x 10" 3 ~ 10" 2 


4.5 x 10~ 3 


In (sin 2 20i 2 ) 


-5.41 ±0.80 


-5.40 


sin 2 26*23 


> 0.8 


0.90 


sin 2 20i3 




1.4 x 10" 3 


TABLE IV. Experimental values versus fit central values for observables using the inputs of 



Table III. Masses are in GeV and all other quantities are dimensionless. Error bars indicate the 



larger of experimental or 1% theoretical uncertainties, as described in the text. 
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FIGURES 
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